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ABSTRACT 

Aims. We present preliminary results of the first Herschel spectroscopic observations of NGC7129 F1RS2, an intermediate mass star-forming 
region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. 

Methods. The PACS instrument was used in line spectroscopy mode (R= 1000-5000) with 15 spectral bands between 63 and 185 pm. This 
provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. 

Results. Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, 
typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. 
Conclusions. It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot 
enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this 
high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of 
UV photon heating and/or non-dissociative shocks. 
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1. Introduction 

, The common paradigm for the formation of low mass stars in- 
■ eludes the spherical infall of a cloud core, the formation of a 
' disk with a protostar at its centre, and a bipolar outflow emitted 
from the inner region of this disk/protostar structure. Large (R 
~ 10,000 AU), low temperature (T ~ 20 K) envelopes, which 
\ are often approximated to be nearly spherical (Terebey, Shu, and 
. Cassen I19841 l are observed around young low mass stars and 
' dominate the observed sub-millimeter emission (Shirley et al 
i2000). Recent observations show additional structures, possibly 
i hot outflow cavity walls heated by UV photons (van Kempen et 
' al [2009 ) and shocks (Gianninni et al |1999i van den Ancker et 
al[20 00 Nisini et al l2002l Arce et al l2007l van Kempen et al 
12010 ). Our understanding of the formation of high mass stars is 
not so mature, although it is clear that there are some signifi- 
cant differences between the low mass and high mass formation 
mechanisms, including the fact that massive stars are predomi- 
nantly born in clusters where interactions between stars are im- 
portant in the formation mechanism. Low mass stars likely form 
in clusters too, but their early evolution may be unaffected by 
neaerby stars until when the first high mass star forms in their 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA 



vicinity. Then the formation process for the low mass stars prob- 
ably changes dramatically and perhaps even ceases. The many 
unanswered questions about these mechanisms include: is there 
a continuum of properties between the low and high mass young 
stellar object; where is the boundary between low mass and high 
mass stars; do these "intermediate mass (IM)" stars form singly 
or in clusters? 

The scarcity of IM stars has made answering such questions 
difficult. High mass stars are also rare in the vicinity of the Sun 
but their high luminosities make it possible to detect and study 
them at much larger distances. Intermediate mass protostars are 
usually defined as those Young Stellar Objects (YSOs) with 
bolometric luminosities between 75 L© and 2 x lO-' L©, although 
both the lower and upper limits of this range vary slightly from 
one author to another There are no examples of well-studied 
very young protostars (e.g. Stage 0/Class 0) with luminosities 
within this range within 400 parsecs of the Sun. There are a 
handful of such objects with distances between 400 and 1500 
pc, which are individually identified and not confused in larger 
OB associations. 

The recent launch of the ESA Herschel Space Observatory 
(Pilbratt et al 2010) provides us with an opportunity to study IM 
stars with previously unavailable molecular probes. The Water In 
Star forming regions with Herschel (WISH) Key Program uses 
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emission from water and other key molecules to probe the con- 
ditions in a variety of star forming regions. 

At a distance of 1260 ± 50 parsecs (Shevchenko & Yakubov, 
[19891 1. NGC 7129 FIRS2 is the most distant of the objects in the 
WISH intermediate mass sub-program. FIRS2 contains a 5 
protostar and has a luminosity of 500 Lq and a hot core (Fuente 
et al l2001ll200 5a). As a very young YSO it was expected that the 
structures dominant in this object would be those of the earliest, 
infalling, stage of the formation process. There is no evidence 
for a large, well-developed disk as one would expect to find in 
an older YSO. On the other hand, a well-developed quadrupolar 
outflow (due to the superposition of two bipolar outflows), has 
been found very close to the protostar (Fuente et al, 2001) . indi- 
cating that some evolution past the collapsing spherical envelope 
stage has occurred. 

A recent analysis of a number of IM YSOs included 
NGC7129 FIRS2 (Grimier et al 2010). This analysis models the 
objects as spherical envelopes, includes all available far infrared 
(Spitzer) and submillimeter (JGMT and IRAM) brightness mea- 
surements. The best-fit model found by Grimier et al (1201 Oi l has 
an envelope mass of 50 Mq, an optical depth at 100 fim of 2.3, 
an inner radius of 100 AU and an outer radius of 18,600 AU. The 
temperature at the inner envelope radius is 300 K and this falls 
to 100 K at a radius of 373 AU where the H2 density is 4.4 x 10^ 
cm""*. The density varies as a power-law with index 1.4. 

We use this model as a starting point for understanding the 
observations described here. We present a first look at NGG7129 
FIRS2, the you ngest k nown IM YSO, with the PAGS instrument 
(Poglitsch et al 120101) on the Herschel Space Observatory. The 
high excitation lines accessible with the instrument should pro- 
vide the best probe available of the warmest, innermost regions. 



2. Observations 

The PAGS instrument was used to obtain data for this project 
during the Herschel Science Demonstration Phase on October 
26, 2009. The pointed line-scan mode was used to obtain data in 
two scans (Obs ID 1342186321 and 1342186322). Data were 
reduced using HIPE v2.4 and then analyzed using IDL. The 
relative and absolute spectral response functions are based on 
ground calibrations and corrected for in-flight observations by 
dividing fluxes below and above 100 pva by 1.3 and 1.1, respec- 
tively Q The continuum levels were found to vary from band to 
band but only a simple polynomial function was required to re- 
move this baseline emission. The noise levels also varied from 
band to band, but a typical value was 0.3 Jy with 0.01 //m chan- 
nel widths in the bands between 100 im\ and 180 /im indepen- 
dent of scan or receiver, which corresponds well to the noise 
levels predicted in HSPOT, the Herschel Science Planning and 
Observing Tool. The spectral resolution is 3,000 at 63.2 /im but 
varies from 1000 to 4000 over the wavelength range used here.. 
The PAGS spectra are imaged in 5x5 pixels, with each pixel 
subtending 9.4" by 9.4" , with no gaps between the pixels. For 
wavelengths shortward of 1 10 ^va this (9.4" ) is the effective res- 
olution but the psf is significantly larger than the pixels at longer 
wavelengths such that at 200 /vm only 40% of the light of a well- 
centred point source falls on the central pixel. The central pixel 
was positioned at the coordinates of the central interferometer 
continuum peak at 2lH3'"1.7' 66°03'23" (J2000). While there 
is still an uncertainty in the flux calibration of approximately 
50%, this does not have a significant impact on the first analysis 
of this object as presented here. 




144.6 1 44.7 1 44.8 1 44.9 1 45.C 
Wavelength [micron] 



179.5 1 80.0 1 80.5 
Wavelength [micron] 



From internal PACS-ICC results. 



Fig. 1. A sample of the spectra detected in these observations 
with PAGS in NGG7129 FIRS 2. 



Table 1 lists 27 lines that were detected, including two pairs 
of blended lines. The fluxes given are measured only in the cen- 
tral pixel. The list of detections include 9 water lines (plus one 
blended with a GO line), 9 unblended GO lines, 3 OH lines (plus 
one blended with a GO line), and 2 [OI] lines. The rest wave- 
lengths, upper level energies as well as the integrated line fluxes 
measured in the central PAGS pixel are provided in Table 1 . A 
number of spectra from these observations are shown in Figure 
1. Note that the line width is essentially all due to the instru- 
ment spectral resolution. Five of these water lines were covered 
with ISO LWS but only upper limits were detected (Grimier et 
al 1201 Oi l. In all five cases the LWS limits are a factor of ten or 
more greater than the detections listed here, with typical limits 
of 10-15 Wm-2. 

In addition to the lines listed in Table 1 there is also a detec- 
tion of the [GII] line at 158 //m and numerous faint lines that can 
not be measured well until the calibration of PAGS improves. 
[G II] 158 pva emission is detected across the entire PAGS field- 
of-view. The other emission lines in Table 1 come from a com- 
pact region, consistent with being smaller than the resolution of 
the instrument. This [GII] emission is likely not directly related 
to the YSO but instead arises from the much larger NGG 7129 
nebula itself and will not be discussed further in this paper. 

3. Discussion 

The GO lines detected in these observations are all lines char- 
acteristic of very hot regions. GO J - 33-32 for instance, arises 
from a level at an energy equivalent to 3093 K, and is not likely 
to be strong unless the excitation temperature is in excess of 
1000 K. The typical excitation temperatures for the other species 
detected here are not as high and so, for this paper, we focus on 
the GO lines and only briefly discuss the other observed lines. 
We first examine our original idea for the dominant structure 
in NGG7129 FIRS2: that this YSO is very young and therefore 
is dominated by a spherical envelope. Is the gas, as traced by 
the most readily observed molecule, GO, compatible with the 
hypothesis of an infalling spherical envelope that can be fitted 
to the dust through SED modelling? We show that this idea is 
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Table 1. NGC 7129 fluxes measured fromPACS spectra 



Species 


Transition 


A 


E„ 


Int. Flux 

10-"* 






fim 


K 




H2O 


O2i2-loi 


179.5 


114.1 


241 




22|-2i2 


180.5 


194.1 


120 




022I-110 


108.0 


194.1 


295 




3o3-2l2 


174.6 


196.8 


227 




P 3l3-2o2 


138.5 


204.7 


194 




p4o4-3i3 


125.4 


319.5 


113 




O4i4-3o3 


113.5 


323.5 


596" 




P 322-2ll 


90.0 


296.8 


113 




o423-3l2 


78.7 


432.2 


308* 




8ig-7o7 


63.3 


1070.7 


66 


OH 


1/2-3/2 


79.1 


181 


209 




1/2-3/2 


79.2 


181 


182 




3/2,7/2-3/2,5/2 


84.4 


291 


290'' 




3/2,7/2-3/2,5/2 


84.6 


291 


181 


01 


'P, - 'P, 


63.2 


227.7 


802 




'Po - 'P2 


145.5 


326.6 


50"^ 


CO 


14-13 


186.0 


580.5 


368 




16-15 


162.8 


751.7 


336 




18-17 


144.8 


945.0 


308 




22-21 


118.6 


1397.4 


273 




24-23 


108.8 


1656.5 


173 




29-28 


90.2 


2400.0 


153 




30-29 


87.2 


2564.8 


246 




32-31 


81.8 


2911.2 


144 




33-32 


79.4 


3092.5 


112 



blended with CO (23-22) Edge effect may have increased the flux 
blended with CO (31-30) very uncertain 



untenable in the face of the data obtained and described here. 
We are led to examine non-spherical geometries to explain these 
observations. There we examine details of the physical and ra- 
diative structure around intermediate mass protostars that could 
prove valuable in understanding star formation. 

3. 1 . CO from a spherical protostellar envelope 

The envelope model used in the analysis by Grimier et al ( 1201 Oi l 
provides a reasonable starting point for comparison to these new 
observations. The envelope in their model extends to 18,600 
AU or 14" at the distance of this object, close to the beam- 
size of these observations. Grimier et al ( 120101 ) do not make 
predictions for the GO lines and so we have reproduced their 
model envelopes and used the RATRAN radiative transfer code 
(Hogerheijde & van der Tak I2000l l to estimate the strength of 
the GO line emission. The model results are plotted in Figure 2 
as a red line showing the integrated line strength versus the up- 
per level rotation The GO line emission observed by PAGS is 
shown on the same plot, indicated by data points with error bars 
showing the 50% uncertainties in the measured fluxes. 

Also shown in Figure 2 are three data points from low-J 
GO, integrated over a 20 arcsec wide aperture. Two GO J-2- 
1 data points are given (from Fuente, private communication), 
the lower one for the velocity range of -1 1 to -9 km/s that may 
be appropriate for the envelope, and the upper one for the entire 
range over which emission is detected (-20 to km/s) appro- 
priate to the outflow. Although we present the -11 to -9 km/s 
ambient cloud velocity value it should be noted that no peak is 
seen in the map in this velocity range. The GO J -3-2 data point 
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Fig. 2. Observed integrated GO line strengths (data points) and 
GO emission calculated from the protostellar envelope model 
presented by Grimier et al (red line), from a GO slab model (blue 
line) as described in text, and the total of the flux predicted from 
these two models (black line). The three points at the left (two 2- 
1 GO and one 3-2 GO) are taken from ground-based observations 
as described in the text. 



is integrated over velocities from -14 to -6.5 km/s (Miskolczi 
et al, I200n and most of the emission is from the outflow - in- 
deed it is not clear that any envelope emission is present in the 
maps shown. A measurement of GO J=l-0 was also made from 
published interferometry images but the flux was too low to ap- 
pear on this plot. Single dish observations of this line show no 
source at this position, as with the GO J=2-l. It is not surprising 
that this low-J GO does not fit the spherical envelope as it likely 
arises from the outflow. 

The structure in the model fit around - 10 is due to two 
effects: the lines are optically thick below that point; and also we 
use a constant beam size of 20" for < 14 and use the actual 
Herschel beam size at higher J. 

The predicted GO emission from the envelope for 7u > 13 is 
many orders of magnitude below the observed emission. This is 
not surprising as the volume within the protostellar envelope that 
exceeds 100 K is very small, with a radius of less than 400 AU. 
This radius corresponds to approximately 0.3" at the distance of 
NGG7129. With a typical Herschel beamsize of 10" for the GO 
lines, this hot region fills a small fraction (less than 10"^) of the 
beam area. We explored other possible spherical envelope mod- 
els using similar constraints to the Grimier et al (1201 Oi l model. 
None of these alternate models have successfully produced a sig- 
nificant amount of GO emission in 7u > 15 lines. 

3.2. CO from other components within the YSO 

A quick analysis to gain some insight into the origin of the high- 
J GO lines is to make a Boltzmann plot (shown as Figure 3 in 
on-line materials). A single temperature LTE component (best fit 
at 52 IK) did not fit well but two components, at temperatures of 
338 K and 799 K, gave a satisfactory fit. However temperatures 
up to 1300 K were found to be consistent within the uncertainties 
in the data. 

One possible source of the high-7 GO lines is a hot core. 
However, the size of such a hot core at the distance of N7129 
is small. This hot core is equivalent to the hot inner part of the 
envelope described above: the size of any such core is very small 
(0.72" X 0.52") and has an upper temperature limit of 890 K 
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(Fuente et al l2005al l. and thus the CO emission will be much too 
small to account for this emission. 

A cavity within the envelope, presumably carved by the ob- 
served outflow, heated by UV photons could present a larger hot 
(T ~ a few 100 K) area within the beam of the PACS observa- 
tions (Spaans et al. 1995 119951 van Kempen et al. 120091 ) Burton 
et al II 9901 show that strong UV emission in clumpy photodisso- 
ciation models can produce significant high-J CO emission but 
in virtually all of their models the mid-J CO lines are several 
orders of magnitude stronger than the high-J CO lines, unlike 
what we see in this object. However one of their models, the 
low UV field, high density, "high temperature solution" (there 
are two possible solutions) can produce high-J lines nearly as 
strong as the mid-J lines. The Burton et al models also predict 
very much stronger emission in lines near J=10-9, providing a 
test of whether or not this model applies to this object. 

While UV heated cavity walls can likely produce the mid-7 
CO line emission that we see (see the similar case for HH46, van 
Kempen et al 2010), the highest-7 CO lines probably require still 
hotter (T > 1,000 K) gas and we suggest that this gas is heated by 
non-dissociative shocks. Here we use the more general RADEX 
code (van der Tak et al 120071) to model CO emitting slabs at the 
temperatures and densities which other studies have shown are 
appropriate for these mechanisms. A detailed analysis to explain 
the CO that considers many other issues will be carried out in a 
future paper. 

Figure 2 shows the results of one such model (blue line). 
We present this as an illustration of how this CO emission can 
arise and stress this is not a fit to the data. Plotted in Figure 2 
are the results from a RADEX model of the CO emission from a 
plane parallel slab with a kinetic temperature of 1 100 K, an H2 
density of 1.0 x 10^ cm""*, and a CO column density of 1.0 x 10'"^ 
cm"^. This temperature is similar to the one derived towards low- 
mass protostars by Neufeld et al 2009 The line width (10 km/s) 
does not affect the result at these low optical depths (typically 
T = 2 X lO"*" in these CO lines). Note that, as with the envelope 
model, the emission is integrated over a beamsize of 15" for 7u < 
14; a beam area appropriate to the diffraction limited Herschel 
beam, is used up to - 20. However for the higher J lines (that 
do not appear in the envelope model) the PACS beam of 9.4" 
is used. The mass contained within the beam, assuming a CO 
abundance of 10""^, is negligible (of order 10"^ M©). 

This slab model provides a reasonable and simple first it- 
eration fit to the observed CO fluxes. The highest-7 CO lines 
provide an excellent match between the observed and predicted 
fluxes. The observed fluxes of the mid-y lines are much greater 
than the modeled fluxes, by a factor of eight for the first of these 
lines (7=14-13). We have not found a single CO component that 
can match the observed flux from these three lines without si- 
multaneously producing significantly too much flux in highest 
lines. It appears likely that at least two separate components are 
required within the current plane -parallel slab geometry in order 
to match all of the observed lines. With further calibration of the 
PACS instrument, the uncertainties in the observed line strengths 
are expected to diminish. As well, this source will be observed 
with HIFI as part of the WISH project. If the mid-7 CO lines 
observed are from shocks then HIFI observations should resolve 
the lines and show them to be broad. A more detailed consid- 
eration of the structure within the protostellar envelope will be 
best undertaken once these additional data and calibrations are 
available. 

We defer a detailed discussion of the other molecules de- 
tected to future papers. The H2O data listed in this paper will be 



combined with recently obtained HIFI data in one of these future 
papers. 

4. Conclusions 

The principal result of these new observations by PACS is that 
the unexpected strength of the observed high excitation (high-7) 
CO lines requires high temperatures as well as relatively high 
densities to be present in the immediate vicinity of NGC 7129 
FIRS2. We have found a single component plane parallel hot 
slab can produce a reasonable fit to the highest-7 CO observa- 
tions but a more complex solution is necessary to fit all the lines 
simultaneously, especially the mid-7 lines. This hot component 
may represent a UV photon heated inner wall of the outflow cav- 
ity or may be related to shocks associated with the outflow itself. 
In the future it may be possible to model all the CO emission us- 
ing a self-consistent 2-D geometric shock and UV heating model 
as has been used in the interpretation of the observations of the 
low mass protostar HH46 (van Kempen et al 1201 Oi l, where sim- 
ilar highly excited CO emission is seen. HIFI observations will 
provide the higher spectral resolution needed to trace the kine- 
matics of these structures. 
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Fig. 3. A Boltzmann plot of the CO emission lines observed with 
PACS. This plot shows the column density versus the excitation 
energy of the upper level. The best fit to a single component 
gives a temperature of 521 K but the fit is not good. Two com- 
ponents with temperatures of 338 K and 799 K produces a better 
fit, as shown. 



